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ABSTRACT 

A geometrically thick dusty torus in NGC 1068 has been unambiguously resolved by an infrared interferometry 
telescope. This implies clouds composing the dusty torus are undergoing supersonic collisions with each other. We 
show that the collisions form strong non-relativistic shocks, which accelerate populations of relativistic electrons. 
Torus reprocesses emission from accretion disk into infrared band. We show that the energy density of the infrared 
photons inside the torus is much higher than that of the magnetic field in the clouds and the seed photons of inverse 
Compton scattering are mainly from the infrared ones. The maximum energy of the relativistic electrons can reach 
a Lorentz factor of 10 5 . We calculate the spectrum of the synchrotron and inverse Compton scattering radiation 
from the electrons in the torus. The relativistic electrons in the torus radiate non-thermal emission from radio 
to 7-ray, which isotropically diffuses in the region of the torus. We find the most prominent character is a peak 
at ~ 0.5- 1 GeV We apply this model to NGC 1068 and find that the observed radio emission from the core 
component S\ can be explained by the synchrotron emission from the relativistic electrons. We predict that there 
is a 7-ray emission with a luminosity of 10 4l) erg/s peaking at ~ 1 GeV from the torus, which could be detected by 
Gamma-Ray Large-Array Space Telescope in the future. This will provide a new clue to understand the physics in 
the torus. The non-thermal radiation from dusty torus may explain the radio emission from Seyfert galaxies. The 
cosmological implications of the non-thermal emission to the 7-ray background radiation have been discussed. 

Subject headings: galaxies: individual: NGC 1068; galaxies: nuclei 



1. INTRODUCTION 

A geometrically thick dusty torus has been recently resolved 
in the well-known Seyfert 2 galaxy NGC 1068 by an infrared 
interferometry telescope (Jaffe et al. 2004; hereafter). The 
dusty torus, as an essential ingredient in the unification scheme 
of active galactic nuclei (AGNs), plays a key role in obscur- 
ing the broad line region (Antonucci & Miller 1985; Antonucci 
1993). How to maintain the thickness of the torus has long 
been puzzling in AGN physics. Krolik & Begelman (1988) 
suggested that the dusty torus is composed of discrete com- 
pressed clouds and supported by the random supersonic mo- 
tion of clouds in z-direction. However, the inevitable collisions 
among the clouds dissipate their kinetic energy, an extra supply 
of kinetic energy to the clouds, such as a star cluster, is thus 
needed (Krolik & Begelman 1988; Jaffe et al. 2004). What are 
the observational consequences of the collisions? 

VLBA (Very Long Baseline Array) has resolved the radio 
components in the core region of NGC 1068. The radio com- 
ponent S\ at the center of the nucleus is coincident with the 
most powerful infrared source in NGC 1068 (Bock et al. 2000; 
Jaffe et al. 2004). The component S\ should therefore relate 
to the dusty torus somehow physically. The inner edge of the 
hypothesized torus will be ionized by UV and X-rays from the 
accretion disk and is expected to be optically thick to free-free 
absorption at GHz frequencies (Neufeld et al. 1994). The free- 
free emission is suggested for the radio emission of the com- 
ponent S\ from hot ionized gas of a temperature T = 6 x 10 6 K 
and number density n = 6 x 10 5 cm" 3 evaporated from dust, and 
the inverted spectrum below 5GHz in NGC 1068 is explained 
by the ionized gas via free-free absorption with a total amount 
of 3400M Q (Gallimore et al. 2004). For typical values of 
the parameters, the ablation rate is of M a bi « 0.49M Q yr" 1 via 
photo-ionization (Krolik & Begelman 1988), the timescale for 



the required mass of the ionized gas will be of f a bi ~ 7 x 10 3 yr. 
However, the free-free emission from hot ionized gas is propor- 
tional to the square of particle density and its cooling timescale 
is of, fft « 3.2T^ 2 n^yT, for the inferred ionized gas, where 
T 6 = T/10 6 K and n 6 = n/10 6 cm" 3 . It is clear that f ff < f aW , 
namely, most of the ionized gas (~ 34OOM ) can not be main- 
tained and will re-condense to dust. So there is no enough ion- 
ized gas to absorb the radio emission below 5 GHz in the model 
suggested by Gallimore et al. (2004). The fact that the Com- 
ponent Si nicely overlaps with the region of the most infrared 
bright torus (Jaffe et al. 2004) strongly implies that the radio 
emission from Si is linked with the torus itself. The explanation 
of this radio emission remains open. Does it relate to collisions 
among the clouds? 

In this Letter, we focus on the observational consequences 
of the collisions among the clouds and show that shocks due 
to collisions accelerate populations of relativistic electrons. 
Non-thermal emission will be radiated as observational con- 
sequences from radio to GeV 7-rays. We apply this model 
to NGC 1068 and suggest that the radio emission of the com- 
ponent Si originates from the synchrotron emission of popula- 
tion of the relativistic electrons. The predicted 7-ray emission 
can be detected by Gamma-ray Large-Array Space Telescope 
(GLAST). This provides a new clue to test the mechanism sup- 
porting the torus. 

2. NON-THERMAL EMISSION FROM TORUS 

2.1. Cloud Collisions and acceleration of electrons 

The detail micro-physics of clouds has been discussed by 
Krolik & Begelman (1988). The clouds are supported by self- 
gravity and magnetic field. Random motions of the clouds 
in the vertical direction support the thickness, H, w Av/tt^, 
where H, is the height of the torus, Av is the random veloc- 
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ity in vertical direction and f^K is the Keplerian velocity of the 
clouds at the mid-plane. We then have H t /R, w Av/v OI b, im- 
plying Aw ~ w or b. It is expected that the clouds in torus are 
undergoing collisions in such a torus. The clumpy torus can 
be described by the dumpiness C = £ c i/ (£ c i)> where E c i and 
(S c i) are the column density of each cloud and the mean col- 
umn density of the torus, respectively. Genzel et al. (1985) 
estimate C w 0. 1 in the Galactic center. The total kinetic energy 
dissipated by the collisions for typical values of the parameters 
is given by Krolik & Begelman (1988) 

£ diss = 1.9 x 10 4 X C E 25 ( - (ergs" 1 ), (1) 



V200km/s 



where the column density of the torus £25 = (£ c i)/10 cm -2 
and the distance to the black hole R pc = R,/lpc. Part of this 
power is used to evaporate dust and left forms shocks which ac- 
celerate populations of relativistic electrons. We do not tackle 
the energy supply to the clouds in this paper. 

The magnetic field in a cloud is expressed in term of the 
plasma parameter defined by (3 = P gas /Pmag, 



- 5 rV2 C V2 s l/2 T3 l/2 w -l /2/? -y2 (Gauss)) 



(2) 
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where TL = H,/R„ f3 . 2 = (3/0.2, C-\ = C/0. 1 and the temperature 
of the torus r 300 = T,/300K. 

The Lamor gyration radius R^ of an electron with Lorentz 
factor 7e is R L = 5.5 x 10~ 8 (75/5-3) pc, where 75 = 7 e /10 5 and 
Z?_3 = B/10~ 3 Gauss. This radius is much shorter than the size 
of the dusty torus. The temperature of the shocked gas during 
the collision is given by 

T S = M 2 T, « 1.2x 10 6 Au 2 (K), (3) 

where March number M. = Av / v^, Av2 = Av /100 km s _1 and 
v± = (kT,/m p ) l l 2 the sound speed. The shock velocity with re- 
spect to the unshocked gas v S h are given by 



w S h = (r+i) 



kT s 



1/2 



J 2.97 xl0 7 r ft 1/2 (cms" 1 ), (4) 

(T-l)m p \ 

where T = 5/3 is the adiabatic index, k the Boltzman constant, 
m p the proton mass and = r s /10 6 K the temperature of the 
shocked gas. The timescale of accelerating electrons reads 
(Blandford & Eichler 1987) 

face = ^-5.8x lO^Zr 1 ^" 1 (sec). (5) 

where c is the light speed. 

The energy density of infrared photons inside the torus, 
which is reprocessed radiation from the accretion disk, can be 
estimated by 



m ir - 



1 .8 x 1 0-^44 W ( er g cm X (6) 
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where L44 = LiR/10 44 erg s _1 is the infrared luminosity from the 
torus. We find m ir » u B , where u B is energy density of the 
magnetic field. Thus the cooling of the relativistic electrons is 
mainly due to inverse Compton scattering of the infrared pho- 
tons. The maximum energy of the electrons will be determined 
by the balances between the shock acceleration and the inverse 
Compton scattering of the IR photons inside the dusty torus. 
The timescale of energy loss due to inverse Compton (IC) scat- 
tering is given by 

3m e c 



4o- T 7 e MiR 



1.7xl0 8 7 J^4l^ c (sec). (7) 



The maximum energy reaches when the acceleration timescale 
is equal to the IC timescale, 

7max w 1.8 x 10 5 L 4 l 4 /2 H l/2 R pc B%X /2 - (8) 
The shock acceleration is known to generally accelerate a 
power law distribution of relativistic electrons as % a , namely, 
the electron distribution function is 



m e c z c ui R 



(9) 
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which follows from the energy equation that the relativistic 
electron is getting energy from the dissipation of the cloud col- 
lision 



7icm e c N(~/ e )dj e = £ e £diss- 



(10) 



where 71c = Co«iR7e (cq = 3.2 x 10 _8 erg _1 s _1 cm -3 ) is the energy 
loss rate of the electron due to inverse Compton scattering (Ry- 
bicki & Lightman 1979). Here £ e is the fraction of the thermal 
energy converted into the energy of the relativistic electrons. 

2.2. Synchrotron Radiation and Inverse Compton Scattering 

For a steady torus, the cloud collisions provide a constant 
energy supply to relativistic electrons. We have shown syn- 
chrotron losses of the electrons are much less than the inverse 
Compton scattering of the infrared photons. However the syn- 
chrotron radiation from the electrons is detectable in radio band 
as one observational consequence. The emergent spectrum 
from the torus will be given by the synchrotron and inverse 
Compton scattering emission. The peak frequency of the syn- 
chrotron emission is 

v syn » 4.2 x 10 13 B_37 5 2 (Hz), (11) 
peaking at middle or near infrared band. This just overlaps with 
the peak frequency of the reprocessing of the disk and thus is 
under-detectable. We approximate the reprocessing of dusty 
torus as a black body with a single black body temperature 

of r eff = (L m /2ir 2 aR 2 Hc) 1/4 w 223L 4 { 4 H^^R^K, where the 
black body radiation constant a = 7.56 x 10 _I5 erg cm" 3 deg" 4 . 
The averaged energy of the infrared photons is ei R w 2.7 'kT e g w 

5.2 x lO^L^TT^flpe^eV in torus, we have the peak of the 
inverse Compton scattering 

£7 » 7i x eiR » 0.52 7? (qr/5.2 x 10" 2 eV) (GeV). (12) 
We use the standard formulations of synchrotron and inverse 
Compton scattering including Klein-Nishina effects to calcu- 
late the spectra from torus for the typical values (Blumenthal & 
Gould 1970). The ratio of energy densities of magnetic field to 
synchrotron photons is of Mfi/wph ~ 10 2 ~ 3 for the typical val- 
ues of the parameters, so self-Compton scattering emission can 
be neglected. Figure 1 shows the spectra for different electron 
index, but the injected energy £ e £diss is fixed. 

The synchrotron spectrum beyond radio band overlaps with 
the reprocessed emission from the torus itself. Since it is much 
lower than the later, this predicted component will be invisi- 
ble or undistinguished. However, the component of > 100 keV 
spectrum will be observable because the overlapped continuum 
from disk/corona has a cutoff at lOOkeV. This is a key feature of 
the present model. This component if detected only originates 
from the torus. 

We emphasize that the non-thermal emission from torus is 
isotropic in the present model. The predicted spectra should be 
testable for any Seyfert galaxies as long as the dusty torus exists 
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FIG. 1 . — The spectrum of component Si . The parameter (3 = 0.2 is suggested in Krolik & Begelman (1988). We use the typical values of B = 5 X 10 3 Gauss, 
7max = HP 

> 7min — 900, «ir — 10 5 erg cm 3 and £ e £diss — 3.0 X 10 40 erg/s. Spectral index a is taken for 2.5 (solid), 2.1 (dashed) and 1.9 (dash-dotted), respectively. 



There is a clear GeV bump predicted from the present model, which reaches a level detected by GLAST. The filled circles, squares and triangles are from Muxlow 
et al. (1996), Roy et al. (1998) and Gallimore et al. (2004). We note that the non-thermal emission from the torus is lower by a factor of 4 ~ 5 orders than the 
emission from the central engine. 



like in NGC 1068. In the present model, the soft X-ray emis- 
sion lines are invisible since they are too faint compared with 
the soft X-ray spectrum of the nucleus. The model of Gallimore 
et al. (2004) predicts strong features in the soft X-ray band, but 
obscured by the torus itself in NGC 1068. It would be inter- 
esting to test the model suggested by Gallimore et al. (2004) 
in Seyfert 1 galaxies, which have dusty tori and the predicted 
luminous soft X-ray emission lines unobscured by the tori will 
be viewed by Chandra. The present model does not depend on 
the types of Seyfert galaxies provided there are geometrically 
thick tori in these objects. So it is not difficult to distinguish the 
two different models in principle. 

3. APPLICATION TO NGC 1068 

The amount of the ionized gas can be estimated from the 
balance between the free-free cooling and the ablation of dust 
in the torus. The ablation timescale is f a bi = M gas /M a bi, where 
the mass of the ionized gas is M gas w l^n^nipT^R] and 
n t h the number density of the ionized gas. We then have 

nth « 11 x 10 4 r 6 1/4 ?r% c 3/2 (cm" 3 ) from t M = t s . The op- 
tical depth of free-free absorption of the ionized gas is ts = 

1.4 x 10" 3 r 6 " 3/2 i^ z «|fl pc , where i/ Gm = j//10 9 Hz and n 4 = 
«th/10 4 cm" 3 . It is thus impossible for free-free absorption 
to cause an invert spectrum at ~ GHz. The free-free emis- 
sion from the ionized gas is of 5.3 x lO^n^T^^z/QHzW 2 /? 3 *. 
(erg s" 1 ), much lower than synchrotron radiation from the rel- 
ativistic electrons. Here the Gaunt factor is neglected in our 
estimations. The thermal electrons from the evaporation has a 
Thomson scattering depth of r es = 2.0 x 10~ 2 n4/? pc . This value 
is consistent with the fraction of the polarized emission to the 
total due to the thermal electrons (Antounucci & Miller 1985). 
Additionally the free-free emission gets a peak flux of 10 38 erg/s 
at <~ O.lkeV, which is overwhelmed by the radiation from the 



core. We can not detect such a faint radiation component. 

The Eddington ratio is roughly L Bo \/L Edd ~ 0.44, where the 
mass of the black hole Mbh = 1O 7 23 M estimated from maser 
observation (Greenhill et al. 1997) and the bolometric luminos- 
ity LboI = 10 44,98 erg/s by integrating the multi wavelength con- 
tinuum (Woo & Urry 2002). It is much higher than the critical 
accretion rate of advection-dominated accretion flow (ADAF). 
This indicates that the radio emission from the Component Si 
can not originate from an optically thin ADAF in NGC 1068. 
Roy et al. (1998) suggested that the component Si may origi- 
nate from synchrotron radiation from the relativistic electrons, 
but it remains open how to accelerate electrons. 

Table 1 Future Observations 

Instrument Energy Band Threshold Note 

(photons cnT 2 s _1 ) 
INTEGRAUSPV 20 keV-8 MeV 2.4 X IF 3 @8MeV x 
INTEGRAL/IB IS 1 20 keV-10 MeV 5.0 X 10~ 5 @ lOMeV X 
GLAST 2 > 100 MeV 1.6 X KT 9 y/ 

Agile 3 30MeV-50GeV 5.0xHr 8 @lGeV ? 

'http://astro.estec.esa.nl/SA-general/Projects/Integral 
2 http://www-glast. stanford.edu/mission. html 
3 http://agile. mi. iasf.cnr.it/Homepage/performances.shtml 

We use the typical values of the parameters in the present 
model. The dissipation rate of the kinetic energy via cloud col- 
lisions is E w 2.0 x 10 41 ergs/s, and £ e = 0.15. The radiation 
luminosity via inverse Compton scattering is of <~ 10 40 ergs/s. 
We find that the radio spectrum of Si can be generally fitted 
by synchrotron emission. The most prominent character in Fig- 
ure 1 is the GeV bump, which is caused by inverse Compton 
scattering of the infrared photons inside the torus. This bump 
corresponds to a flux of ~ 10~ 9 photons cm" 2 s" 1 at 1 GeV and 
~ 10~ 8 photons cm" 2 s" 1 at > lOOMeV. Detection of this com- 
ponent is essential to probe physics inside the torus. 

Seyfert galaxies have been observed by EGRET, showing up- 
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per limit with 2cr (Lin et al. 1993). Table 1 lists the sensitiv- 
ities of several missions. The last column gives the comments 
to observations. We find that the Italian mission Agile could 
marginally detect GeV bump in NGC 1068 and GLAST has ca- 
pability of detecting this GeV bump. 

The synchrotron self-absorption (SSA) depth is numerically 
given by r syn « 1.5 x 10~ 5 K 56 R^ (a = 2.5, 7min = 900, B = 
5.0 x 10" 3 Gauss) at 1 GHz, where K 56 = K/10 56 , for a spatially 
homogeneous distribution of the relativistic electrons inside the 
torus. This shows less importance of the SSA, however we note 
that this may severely underestimate r syn depending on the spa- 
tial density distribution of the electrons. This is because we use 
the number density of the relativistic electrons averaged over 
the entire torus, which is much lower than the local. Addition- 
ally, the GHz spectrum may rely on the break energy of the 
electron distribution, which is determined by the detail of pro- 
cesses of the acceleration and cooling. We remain it as a free 
parameter to fit the spectrum since we pay main attention to 
the observational consequences of cloud collisions. Future de- 
tail numerical simulations of cloud collisions may display some 
details of accelerations so that the minimum energy of the elec- 
trons and SSA depth will be given. 

4. DISCUSSIONS 

Collisions simultaneously lead to loss of angular momen- 
tum of molecular clouds and then determine the accretion rate 
of the black hole (Krolik & Begelman 1988). Therefore the 
present model predicts a strong correlation between the radio 
emission from torus and a parameter representing the accretion 
rate. Indeed Ulvestad & Ho (2001) find a very strong correla- 
tion, Pfc m oc irn (estimated from their Fig 5) in the Palomar 
Seyfert galaxies, where L[o m] is the [O III] luminosity and Pecm 
is the radio emission power at 6cm. L[o mj is a good indicator 
of the accretion rate of the black hole since L[o m] oc L loa , where 
Li on is the ionizing luminosity. We note that the collisions lead 
to angular momentum loss of clouds, which will be accreted 
onto the black hole. Therefore the radio emission and L[oni] are 
related with the physics of cloud collisions. The P6cm-L[oni] 



correlation could be explained by the present model. The pre- 
dicted radio flux in the present model depends on electron index 
a, a more careful study on this will be carried out in future. 

The isotropic radiation of the relativistic electrons may imply 
a significant contribution to the 7-ray background radiation. If 
the typical 7-ray luminosity from torus is of L 7 GN 10 40 erg/s, 
we have the contributed flux to the 7-ray background radiation 
F 7 w A gnL 7 GN AV 1 '4ndl, where c/> A gn is AGN number den- 
sity peak, AV is the shell volume of the AGN number peak 
at redshift z ~ 1 and di is the distance the peak shell. We 
have the contribution of AGN torus to the 7 -ray background 
F 1 w 0.6 keV cm" 2 s" 1 sr" 1 for </> A gn = 5.0 x 10" 5 Mpc" 3 from 
X-ray observation (Steffen et al. 2003) and the peak shell 
width AR w d^. The deduced 7-ray background emission is 
of 1 .0 ke V cm" 2 s" 1 sr" 1 below 1 GeV (see Figure 1 in Loeb 
& Waxman 2000). The predicted GeV flux significantly con- 
tributes to background emission below 1 GeV. It is interesting 
to note that the background emission spectrum shows increases 
below 0.6 GeV (Loeb & Waxman 2000). It is worth investigat- 
ing this subject in detail further. 

5. CONCLUSIONS 

We suggest that the collisions among the clouds inside the 
torus produce populations of relativistic electrons in active 
galactic nuclei. Non-thermal emission from relativistic elec- 
trons is detectable in radio and 7-ray band. The radio emission 
from the Component Si in NGC 1068 can be explained by the 
present model. It is predicted that there is a GeV bump in NGC 
1068, which is caused by inverse Compton scattering off the in- 
frared photons inside dusty torus. This GeV component can be 
detected by GLAST. We also show that the GeV bump of radia- 
tion from torus could significantly contribute to the background 
radiation of < 1 GeV. 

I thank the anonymous referee for the helpful and prompt 
comments. This research is supported by Grant for Distin- 
guished Young Scientists from NSFC, Hundred Talent Program 
of CAS and 973 project. 
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